JC-98-06 



|cond-mat/9812416 



Critical Exponents near a Random Fractal 

Boundary 

John Cardy 

Department of Physics 
Theoretical Physics 
1 Keble Road 
Oxford 0X1 3NP, UK 
& All Souls College, Oxford 



Abstract 

The critical behaviour of correlation functions near a boundary 
is modified from that in the bulk. When the boundary is smooth 
this is known to be characterised by the surface scaling dimension 
x. We consider the case when the boundary is a random fractal, 
specifically a self-avoiding walk or the frontier of a Brownian walk, in 
two dimensions, and show that the boundary scaling behaviour of the 
correlation function is characterised by a set of multifractal bound- 
ary exponents, given exactly by conformal invariance arguments to be 
A n = ^(Vl + 24nx + ll)( v / l + 24nx - 1). This result may be in- 
terpreted in terms of a scale-dependent distribution of opening angles 
a of the fractal boundary: on short distance scales these are sharply 
peaked around a = 7r/3. Similar arguments give the multifractal ex- 
ponents for the case of coupling to a quenched random bulk geometry. 



The subject of boundary critical behaviour is by now well under- 
stood, particularly in two dimensions 0. The two-point correlation function 
(4>(r)(f>(R)) of a scaling operator <f), which behaves in the bulk at large dis- 
tances at the critical point as \r — R\~ 2x , where x is the bulk scaling dimension 
of </>, is modified when one of the points (say r) is close to the boundary to 
the form 

(0(r)0(i?)) ~ \r\- x \R\- x \R/r\-*, (1) 

where x is the corresponding boundary scaling dimension, and the angu- 
lar dependence has been suppressed for clarity. In two dimensions, the role 
played by x is emphasised by making the conformal mapping z — > In z of the 
upper half plane to a strip of width tt: in that geometry the correlation func- 
tion decays exponentially along the strip with an inverse correlation length 
equal to x ||. 

Eq. |l| refers to the case when the boundary is smooth (at least on scales 
<C r) and it is an interesting to ask whether these results are modified when 
the boundary is a fractal on these scales. The example of an edge (or corner 
in two dimensions) on the boundary was analysed some time ago [Q and it 
was shown that new edge scaling dimensions arise which depend continuously 
on the opening angle a. In two dimensions || this dependence is given by 
conformal invariance arguments by the simple form x(a) = irx/a. This sug- 
gests that close to a fractal boundary, which may be thought of as presenting 
a distribution of opening angles (which perhaps also depends on the scale at 
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which it is probed), an even more complicated behaviour should obtain. 

In the case of a random fractal, one also expects to see behaviour char- 
acteristic of correlation functions in a quenched random environment. That 
is, they may exhibit multiscaling, which means that the average of their nth 
power does not scale in the same way as the nth power of their average. 
In this letter, we consider two cases where this problem is exactly solvable 
using conformal invariance methods in two dimensions, namely when the 
fractal boundary is a self-avoiding walk, and when it is the frontier (exte- 
rior boundary) of a Brownian (ordinary) random walk. In fact, both cases 
turn out to give identical results. Our methods are a simple generalisation of 
arguments due to Lawler and Werner 0, who have derived exact relations 
between multifractal exponents corresponding to self-intersection properties 
of Brownian walks in two dimensions. This corresponds to the special case 
when is a free scalar field satisfying Laplace's equation. This physically 
interesting example, and its relation to the exponents of star polymers, was 
in fact discussed some time ago in 4 — e dimensions by Cates and Witten 
[[?]]. The results of Lawler and Werner have recently been given an elegant 
interpretation and derivation by Duplantier || in the context of coupling the 
system to a randomly fluctuating metric. 

Consider for definiteness self-avoiding walks 7 which are constrained to 
pass through the origin O. In order to be able to apply conformal invariance 
arguments, we work in the fixed fugacity ensemble, in which each walk of 
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length L is counted with a weight y L , at the critical point where y~ x = pi, 
the lattice-dependent connective constant. The properties of the measure 
on walks on distance scales much larger than the lattice spacing are then 
supposed to be conformally invariant. Denote the radial coordinate by p. 
We want to focus on those walks which have a typical linear size R, and for 
which the origin O is a typical interior point. Without loss of generality for 
computing scaling dimensions, we may then restrict the walks 7 to have the 
form of a pair of mutually avoiding self-avoiding walks, starting from O and 
ending on the circle p = R. In the same spirit we may take these points to 
be the first intersections of the walks with this circle. The region bounded 
by 7 and an arc of the circle p = R is thus simply connected. In this region 
we consider a critical system (for example an Ising model) with a suitable 
conformally invariant boundary condition on 7 and on p = R (for example, 
that the spins are free). Consider the correlation function (<p(r)<j)(R')), where, 
without loss of generality for computing scaling dimensions |7| we can choose 
|.R'| ~ R, and we are interested in the limit where r <C R- The geometry is 
illustrated in Fig. |I|. 

Obviously this correlation function depends on 7, but we may hope to be 
able to compute suitable averages of this quantity over realisations of 7. By 
analogy with the case of a smooth boundary, we expect that 

((f)(r)<j){R')) n ~ \r\- nx \R\- nx \R/r\- Xn , (2) 

where the overline means an average over all the allowed realisations of 7. 
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Figure 1: Geometry in which the simply connected region RiOR 2 is bounded 
by a self-avoiding walk 7 and an arc of the circle p = R. This region contains 
a critical system of which the correlation function of local operators at points 
r and R' is of interest. This region, excluding the disc p < \r\, is conformally 
mapped into a long strip of width n and length £ e g (7, r/R). 
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Note that if the average were over smooth boundaries only, we would expect 
A ra = nx, with only the prefactor modified by the averaging. 

We now make a conformal transformation which maps the fractal bound- 
ary into a smooth one. It is convenient to exclude the disc p < r. This leaves 
the simply connected region bounded by two segments of 7 and arcs of the 
circles p = r and p = R. By the Riemann theorem, the interior of this re- 
gion may be mapped conformally by an analytic function z — > f(j, \r\, R; z) 
onto the interior of a strip of width ir, but with a length which is not simply 
I ~ \n(R/r), but which will also depend on 7. Let us denote this by £ e s{l, 
The correlation function (<p(r)<f)(R')) will be related by this conformal map- 
ping to one between operators located near the ends of this strip. Taking the 
nth power and averaging we see by comparison with (g) that 

At this point, we need further information about averages of quantities 
of the form e~ p£cB , for arbitrary p. In particular, let us consider the case 
where n — 1, and 0(r) is the M-leg operator for M mutually self-avoiding 
walks. It is convenient to generalise slightly and take <f)(R') to be a product 
of distinct single leg operators corresponding to the walks all ending on the 
arc of the circle p = R. The correlation function then gives the number 
of such walks which all begin at r and end at a distance ~ R, and, in the 
critical fugacity ensemble, scales in the bulk like |r — R\~ x m~ Mx ^ ; and near 
a smooth boundary according to (JTJ) with a boundary scaling dimension Xm- 
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Coulomb gas arguments |9j lead to the conjectures % = j$M 2 — i and 
xm = |M(3M + 2), which have been confirmed by numerical work and 
various other known exact information. If we now imagine taking an M-leg 
star polymer near the fractal boundary 7, which itself is a 2-leg star polymer, 
and performing the same average over the realisations of 7, the result will be 
an (M + 2)-leg star polymer in the bulk. We conclude that 
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-(x M+2 -x 2 )e. ^ e -x M tegM) ; (4) 



where the factor e X21 arises from the normalisation of the probability dis- 
tribution of 7. Note that although this is initially defined only for M a 
non-negative integer, it may be continued to other real values for which the 
average exists. Comparing with (|3|) we may therefore choose M such that 
Xm — nx, solve for M using the exact conjecture for % given above, and 
substitute this into the exact form for xm+2- After some simple algebra, this 
gives the result for A n quoted in the abstract. 

A similar argument may be made when 7 is the exterior boundary of a 
Brownian walk. In this case the measure is rigorously known to be confor- 
mally invariant. As in the example considered by Lawler and Werner |J, 
one may now invoke the exact conjectures made by Duplantier and Kwon 
1(J for the relevant dimensions of M mutually avoiding ordinary random 



walks. However, the final result is identical, indicating that there is strong 
element of universality between these two fractals, not only with respect to 
their fractal dimensions fllTf . 
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The nontrivial dependence of A n shows that correlation functions near the 
boundary have a broad distribution of values. In particular, their average 
value, which scales as (r/R) Xl , may be quite different from a typical value. 
As argued, for example, in Ref . fl2| , the typical dependence should be of the 
form (r/R) x \ where A' = d\ n / dn\ n= o = 3x. Interesting enough this is the 
behaviour which would obtain in wedge of interior angle a = ir/3. This idea 
may be made more explicit by interpreting (|3|) in terms of an average over a 
scale-dependent distribution P(a,£) of interior opening angles: 

/ da P(a, t) e -( n ^ nM ~ e~ Xnl (5) 
Jo 

Requiring that this be valid for all positive real n determines the form of P. 
First we see that the behaviour of A ra ~ nx/2 as n — > oo at fixed I implies 
that the effective angle in this regime is a ~ 2tt. This is in agreement with a 
general argument of Cates and Witten 0. If we set uj = - — | and u = nx, 
and define P(u,£)duj = P(a,£)da the above equation simplifies to 

J™ du P(u, t) e- uul ~ exp (-(5£/24)( v / l + 24m - 1)) . (6) 

Making the ansatz P(u,£) ~ e M//24 e -^( aw + 6 /^) anc i us i n g steepest descent 
then leads to the solution 

P^^expf-l^-^) 2 ) , (7) 

valid for large t (i.e. r/iJ < 1) and where we have suppressed more slowly 
varying prefactors. Note that this result is independent of x, consistent with 
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it being an intrinsic property of the fractalQ. It shows that the effective 
opening angle has a broad distribution which, however, becomes more and 
more strongly peaked around the typical value lu — | (a — |) as r/R — ► 0, 
with a width of order (ln(i?/r)) -1 / 2 . 

A similar calculation may be carried out when the point O is the root of 
an iV-leg star polymer, by replacing xm+2 — %2 in @ by xm+n — xn- The 
case N = 1 gives the end multifractal exponents, which, as first pointed out 
by Cates and Witten, are different form those which arise when O is a typical 
interior point. 

In the case considered by Cates and Witten 0, where is a massless 
scalar field (with bulk scaling dimension x = 0) satisfying Dirichlet conditions 
= on the boundary, the appropriate boundary scaling operator is d±(p 
with dimension x — 1. In that case we find that A(I) = | = 2 — D, where 
D — g is the fractal dimension of the boundary. This is a consequence of 
the fact that <fi satisfies Laplace's equation, equivalent to the conservation 
of particle flux in the Brownian interpretation @. Note that we were lead 
to this unique result from making simple assumptions (rigorously grounded 
in the Brownian case) about the conformal invariance of the measure on 
7. This suggests that all such curves which, with probability one, bound a 

simply connected region when viewed on macroscopic distance scales will fall 

1 However, it should be noted that the method of averaging used here, which sums over 
all realisations of 7 passing through a given point O at given distance r from a fixed point, 
tends to emphasise those parts of the fractal for which the opening angle is small. 
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into this universality class and, in particular, will have D — |. 

One of the most interesting features of our main result is that the X n 
are generally, even for integer n, irrational (but algebraic) numbers. This 
is not in disagreement with any established results, since even if the bulk 
critical theory is unitary, there is no reflection positivity in the presence of 



a fractal boundary and so the theorem of Friedan, Qiu and Shenker |L3 



is evaded. However, most examples of exactly calculable critical exponents 
in two dimensions have, even in nonunitary cases, led to rational values. 
Recently Duplantier || has given an interesting interpretation of these type 
of results in the case of a general mixture of Brownian and self-avoiding 
walks, by considering the effects of coupling the system to a fluctuating 
background metric (quantum gravity). He was able to argue that, just as 
the scaling dimensions of overlapping objects in flat space should be added 
to obtain that of the composite, when they are coupled to quantum gravity 
their dressed scaling dimensions are additive if they avoid each other. In 
this way, by going back and forth between flat space and quantum gravity, 
and using the relation between ordinary and dressed scaling dimensions first 
obtained by Knizhnik, Polyakov and Zamolodchikov (KPZ) ||14|| , specialised 
to the case c = appropriate to Brownian walks, he was able not only to 
recover the results of Lawler and Werner ||, but also to derive the earlier 
conjecture of Duplantier and Kwon |L(J. Thus, for this example, his methods 
are more powerful than the simple arguments we have used above, since we 
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found it necessary to invoke the conjectured values for the M-leg scaling 
dimensions. Since our main result is a simple generalisation to the case when 
x ^ 1, one would expect that similar quantum gravity methods might apply. 
However, our result is supposed to be valid for a bulk theory with arbitrary 
central charge c, and it is therefore not clear why the KPZ relation with c = 
should appear in this more general case. 

We have given an exact formula for the multiscaling boundary exponents 
of an arbitrary conformally invariant two-dimensional critical system close to 
a random fractal boundary. This is the first example when such a multifractal 
spectrum with a non-trivial analytic structure has been found exactly. The 
basic method was to realise that this kind of geometric quenched disorder may 
be gauged away by making a suitable conformal transformation, at the cost 
of modifying the moduli (in this case £ = \n(R/r).) The effective distribution 
of £ e g is then probed by replacing the critical system by one with c = (in 
our case, self-avoiding walks) for which the partition function is unity and 
therefore the quenched average of a correlation function is the same as its 
annealed average, which is more simply dealt with. 

Similar ideas may be applied to a critical system coupled to bulk quenched 
disorder in the form of a random metric (which may be realised as the con- 
tinuum limit of a randomly connected lattice.) In this case we may consider 
an annulus of inner and outer radii r and R respectively, which in the case 
of a flat metric may be mapped conformally by the transformation z — > In z 
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to a flat metric on a cylinder of perimeter 27r and length I = ln(R/r). The 
inverse correlation length along this cylinder is then equal to the bulk scaling 
dimension x ||. An arbitrary metric g on the annulus is also conformally 
equivalent to the cylinder with a flat metric but with a length £ e d(g,£). In 
analogy with @ the multifractal bulk exponents A^ , which govern the decay 
of the quenched average of the nth power of the bulk correlation function, 
are given by 



e - X n £ rv, e -nxe cS (g,e) (g) 

We now consider the special case when the critical system case has c = 
0. In that case the quenched and annealed averages are identical, and the 
respective scaling dimensions X Q and X of an operator in flat space and when 
coupled to a fluctuating metric are related by the c = version of the KPZ 
relation |14| Xq = + X). Thus if we now set Xq = nx and solve for X, 



this will yield A^ . The result is 

\ b n = \ (Vl + I2nx - l) . (9) 



This result for the scaling dimensions when coupled to a quenched random 



lattice, was derived for the case n = 1 by Baillie, Hawick and Johnston [[15 
but in fact these correlation functions exhibit multiscaling, with a whole 
spectrum of such exponents. Note that in this case the typical decay of a 
correlation function is determined by (A 6 )' = 3a;. 

The author thanks G. Lawler and B. Duplantier for explaining their ideas, 
and the Fields Institute, Toronto, where this work was started, for its hospi- 
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After this work was completed preprints by Duplantier []16| and Aizen- 



man, Duplantier and Aharony [|17| appeared in which, among other things, 
it is argued that the fractal dimension of the accessible perimeter of a per- 
colation cluster has fractal dimension D = ~, consistent with the general 
arguments advanced above. 
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